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Our  findings  are  that: 

1.  A more  accurate  initial  position  and  track  approaching  the  accuracy  of 
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3.  The  current  method  of  analyzing  the  initial  stream-function  field  in  the 
area  influenced  by  the  storm  circulation  results  in  a loss  of  useful  observations 
(mainly  satellite-derived  wind  estimates)  in  the  outer  fringes 
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of  the  current  SANBAR  forecast  grid  are  producing  serious  errors  beyond  24 
hours  range  in  some  storms  in,  or  recurving  info,  the  westerlies  of  middle 
latitudes . 

5.  Some  storm  tracks  respond  to  baroclinic  effects  in  the  large-scale 
features  of  the  surrounding  region.  In  some  additional  instances  the  storm 
track  may  be  directly  influenced  by  asymmetric  baroclinic  effects  in  the  storm 
structure  itself. 
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INTRODUCTION 


The  goal  of  this  study  is  to  establish  the  value  of  satellite  observations 
for  the  prediction  of  tropical  cyclones.  We  are  especially  concerned  with 
observations  which  can  be  derived  with  ease  from  the  McIDAS  system  (U.  of 
Wisconsin.  1973)  at  the  Air  Force  Geophysics  Laboratory  (AFGL).  More- 
over, we  are  particularly  concerned  with  prediction  of  the  cyclone  tracks  by 
SANBAR  , (Sanders  et^aL  1975),  a barotropic  dynamical  model  applied  to  the 
large-scale  flow  pattern  surrounding  the  storm,  averaged  through  the  depth 
of  the  tropical  troposphere.  This  model  is  used  operationally  at  the  National 

Hurricane  Center  (NHC). 

Because  of  lack  of  adequate  data  within  the  area  influenced  by  the  storm, 
the  analysis  within  this  region  is  model&d  on  the  basis  of  the  estimated  size, 
intensity,  and  current  displacement  vector  of  the  storm,  as  described  by 
Sanders  et  al.  (1975).  Elsewhere  the  analysis  is  based  on  winds  observed 
by  rawinsondes,  by  aircraft.and  by  ships;  and  (recently)  on  wind  estimates 
based  on  cloud-motion  vectors  derived  from  geosynchronous  satellites. 

Our  work  has  focused  on  improving  the  accuracy  of  the  SANBAR  fore- 
casts through  more  accurate  determination  of  the  position  and  track  of  the 
storm,  and  through  improvement  of  the  large-scale  analysis  by  more  exten- 
sive use  of  the  satellite  wind  estimates.  In  this  report,  we  discuss  the  first 
of  these  efforts  in  detail,  and  examine  intensively  some  particularly  poor 
operational  forecasts  to  determine  the  cause  of  the  errors  and  to  estimate 
the  likelihood  of  improvement  by  increased  use  of  the  satellite  data.  All 
cases  were  taken  from  the  1975  hurricane  season  in  the  Atlantic  Ocean,  the 
Caribbean  Sea. and  the  Gulf  of  Mexico.  The  tracks  of  these  storms,  deter- 
mined by  careful  analysis  after  the  fact  (the  "best  tracks"),  are  shown  in 
Fig.  1.  The  real-time  performance  of  SANBAR,  of  several  other  guidance 
forecasts,  and  of  the  ultimate  "official"  prediction  issued  by  NHC,  is  shown 
in  Table  1.  The  various  other  guidance  forecasts  are  discussed  by  Newmann 
<?t  al.  (1972),  Newmann  and  Hope  (1973),  and  Newmann  and  Lawrence  (1975). 


TABLE  1.  Displacement  errors  (nautical  miles)  for  1975  guid- 


ance forecasts. 
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Figure  1.  Best  tracks  of  1975  storms 
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USF  OF  BEST  TRACK  POSITIONS  AND  DISPLACEMENT  VECTORS 

The  McIDAS  system  with  its  superior  navigation  offers  the  possibility 
of  more  accurate  tracking  of  a tropical  storm  up  to  the  initial  time  of  the 
forecast.  This  determination  could  be  based  on  application  of  the  correla- 
tion technique  (U.  of  Wisconsin,  1973)  either  to  the  storm  cloud  system  as 
a whole,  or,  preferably,  to  the  eye  if  it  is  discernible.  To  simulate  the 
best  of  all  possible  worlds,  we  have  taken  the  best-track  information  for  74 
of  the  78  cases  in  which  operational  SANBAR  forecasts  were  made  in  197  5.1 
These  forecasts  were  rerun  by  NHC  with  the  operational  data  base  but  with 
the  best-track  initial  positions  and  displacement  vectors  instead  of  the  oper- 
ational ones.  The  operational  and  best-track  values  are  given  for  each  case 
in  Table  2.  The  best-track  vector  was  obtained  from  the  6-hour  displace- 
ment immediately  prior  to  initial  time,  to  simulate  the  real-time  situation. 

Position  errors  for  the  operational  forecasts  and  for  the  reruns  are 
summarized  in  Tables  3 and  4,  respectively.  Verification  in  both  instances 
was  taken  from  best-track  positions.  The  number  of  verifications  decreased 
monotonically  with  forecast  range  as  increasing  numbers  of  tracks  passed 
inland  or  eastward  out  of  the  range  of  NHC  forecast  responsibility.  In  a few 
instances  the  storm  weakened  within  the  NHC  area  so  as  to  obviate  the  ne- 
cessity of  tracking.  It  is  unlikely  that  the  accuracy  of  the  best-track  infor- 
mation could  even  be  achieved  in  real  time,  if  only  because  the  position  at 
a given  time  is  determined  by  observed  fixes  after,  as  well  as  before,  that 
time.  In  real  time  we  will  be  hard  put  to  know  whether  the  most  recent 
change  in  track  represents  a significant  change  in  the  displacement  rate  of 
the  storm,  or  merely  a short-period  excursion.  Only  analysis  after  the 
fact  can  effect  that  determination.  Therefore,  the  results  in  Table  4 repre- 
sent a limit  of  optimism  in  this  respect. 

Subject  to  the  above  qualifications,  the  results  show  a substantial  im- 


1 The  four  remaining  forecasts  could  not  be  rerun  for  a variety  of  reasons 
unrelated  to  the  accuracy  of  the  operational  predictions. 
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TABLE  2 Operational  and  best-track  initial  positions  and  displacement  vectors. 
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TABLE 

3 . Me  an 

position  < 

errors  (nautical 

miles)  for 

operational 

forecasts 

00  hr 

12  hr 

24  hr 

36  hr 

48  hr 

72  h r 

Amy 

8 

72 

107 

117 

197 

338 

(7) 

16) 

15) 

15) 

(5) 

13) 

Blanche 

16 

53 

101 

135 

-- 

-- 

(4 ) 

13) 

12) 

(1) 

Carol  ine 

7 

33 

73 

116 

152 

-- 

15) 

(4) 

13) 

12) 

11) 

Doris 

1 1 

72 

104 

174 

240 

369 

110) 

19) 

18) 

17) 

16) 

'4) 

E loise 

18 

60 

95 

134 

198 

28  3 

115) 

114) 

M3) 

(12) 

ill) 

'9) 

Faye 

19 

73 

151 

218 

268 

483 

114) 

(13) 

ill) 

19) 

17) 

1 5 ) 

Gladys 

11 

76 

144 

233 

347 

459 

117) 

(17) 

116) 

115) 

114) 

112) 

Bailie 

46* 

65 

— 

— 

— 

— 

12) 

11) 

Mean 

1 5 

67 

121 

181 

261 

39  3 

174) 

167) 

158) 

(51) 

'44) 

133) 

- Significant  at  10%  level 
0 Number  in  sample 


% 
r i 

TABLE  4. 

Mean  position  errors 

(nautical 

miles) 

for  SANBAR) 

reruns. 

"1 

i 

OOhr 

12hr 

24hr 

36h  r 

48hr 

72hr 

Total 

Amy 

0 

63 

111 

155 

237 

400 

(7  ) 

16) 

15) 

15) 

(5) 

(3) 

(24) 

i 7 ] 

[3] 

[3] 

[2] 

[1] 

[1] 

[10] 

Blanche 

0 

27 

72 

137 

14) 

13) 

12) 

ID 

16) 

[4] 

[3] 

[2] 

[0] 

[5] 

Caroline 

0 

24 

42 

79 

220 

15) 

14) 

(3) 

(2) 

ID 

110) 

[5] 

[2] 

[2] 

[1] 

[0] 

[5] 

Doris 

0 

42 

74 

112 

171 

302 

■ 

110) 

(9) 

(8) 

(7) 

(6) 

(4) 

(34) 

[10] 

[7] 

[6] 

[6] 

[5] 

[4] 

[28] 

Eloise 

0 

41 

68 

105 

179 

310 

115) 

(14) 

113) 

112) 

(ID 

(9) 

159) 

[15] 

[9] 

[10] 

[7] 

[6] 

[5] 

[37] 

[ ' 
i v 

Faye 

0 

67 

136 

181 

195 

411 

114) 

(13) 

(11) 

(9) 

(7) 

(5) 

(45) 

[14] 

[10] 

[6] 

[6] 

[6] 

[4] 

[32] 

Gladys 

0 

52 

122 

201 

293 

430 

117) 

(17) 

(16) 

(15) 

(14) 

(12) 

(74) 

[17] 

[10] 

[10] 

[9] 

[8] 

[6] 

[43] 

i k 

Hal  lie 

0 

88 

_ _ 

r *i 
1 

12) 

(1) 

ID 

[2] 

[0] 

[0] 

f * 
1 » 

Mean 

0 

50 

99 

152 

224 

376 

\ + 

(74) 

(67) 

(58) 

151) 

(44) 

133) 

1253) 

% 

1 ).  * 
. V 

[74] 

[44] 

[39] 

[31] 

[26] 

[20] 

[160] 

: >•» 

Change  from  operational  forecasts 

: % 

-15 

-17 

-22 

-29 

-37 

-17 

* 

Total  excluding  OOhr 

* 

1 ) Number 

in  sample 

I Ik 

[ ] Numbe 

r in  sample  which  had 

position  errors  equal  to  or 

less  than  those 

| r. 

of  operational  forecasts. 
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provement  at  all  ranges.  From  previous  experience  we  believe  that  the  ben- 
efit due  to  correction  of  the  initial  position  alone  would  be  a maximum  at  the 
initial  time  and  would  have  largely  vanished  48  hours  later',  Since  we  find 
here  the  greatest  improvement  in  the  range  from  36  to  48  hours,  it  appears 
that  the  corrected  displacement  vector,  which  immediately  improves  the 
stream -function  analysis  in  the  region  influenced  by  the  storm  at  the  initial 
time,  has  also  a surprisingly  beneficial  large-scale  effect  at  later  times.  By 
72  hours,  however,  the  improvement  is  decreasing  rapidly. 

We  note  that  the  mean  errors  for  any  forecast  range  for  any  individual 
storm  differ  only  insignificantly  from  the  mean  over  all  storms.  That  is, 
none  of  the  1975  storms  were  exceptionally  difficult  (or  easy)  to  predict.  This 
uniformity  is  probably  due  to  the  lack  of  any  notably  eccentric  tracks  during 
this  year,  as  can  be  seen  by  inspection  of  F’ig.  1.  In  the  event  of  loops,  cusps, 
halts,  and  sudden  starts,  it  is  likely  that  the  best-track  information  would 
have  differed  from  the  operational  data  more  widely  than  shown  in  Table  2. 

It  is  thus  probable  that  the  improvement  in  forecast  accuracy  would  be  greater' 
In  this  respect,  then,  the  improvement  shown  in  Table  4 may  understate  the 
potential  of  the  best-track  information. 

Two  aspects  of  the  displacement  error  vectors  were  examined,  for  both 
the  operational  and  the  rerun  forecasts.  The  first  of  these  was  simply  the 
speed  error,  expressed  in  nautical  miles  per  24  hours.  The  second  was  a 
"direction  error",  the  perpendicular  distance  from  the  forecast  position  to 
the  observed  displacement  vector.  The  definitions  of  these  errors  is  illus- 
trated by  the  sketch  in  Fig.  2.  This  was  done  separately  for  the  first  ami 
for  the  second  24-hour  periods  of  each  forecast.  Results  are  summarized 
for  mean  magnitudes  in  Table  5 and  for  algebraic  means  in  Table  6.  In  the 
latter,  a negative  value  means  too  slow  a forecast  speed,  or  a forecast  posi- 
tion to  the  left  of  the  observed  track,  as  shown  in  Fig.  2 

In  Table  5 we  see  that  speed  errors  and  direction  errors  are  of  compar- 
able magnitude,  that  both  grow  substantially  from  the  first  to  the  second  24- 
hour  period  of  the  forecast,  and  that  both  types  are  reduced  slightly  for  both 
periods  in  the  rerun  forecasts.  The  statistical  significance  of  the  imprnve- 
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Speed  error 


Predicted 


Figur 


Sketch  illustrating  spued  and  direction  errors 


TABLE  5.  Mean  magnitudes  of  speed  and  direction  errors. 


Number  in  sample 


ment  in  the  overall  means  is  not  large  , but  inspection  of  Table  5 shows  that 
improvement  occurred  for  the  individual  storm  means  in  about  three-quarters 
of  the  instances. 

In  Table  6,  some  interesting  biases  appear.  The  forecast  speeds  are 
slow  and  the  displacements  are  to  the  left  of  the  observed  tracks.  Surpris- 
ingly, the  biases  are  not  reduced  in  the  rerun  predictions  based  on  best-track 
information.  Reference  to  Fig.  1 shows  that  the  1975  storms  tended  strongly 
to  move  along  clockwise-curved  tracks  (i.e.  to  recurve)  and  accelerate  during 
the  period  when  forecasts  were  made.  This  tendency  was  especially  pro- 
nounced for  Eloise,  Faye,  and  Gladys,  for  which  the  majority  of  predictions 
were  made.  < In  31  of  the  74  cases  the  initial  position  was  north  of  30°N.)  Evi- 
dently the  forecast  tracks  did  not  show  these  characteristics  to  a sufficient 
degree,  presumably  because  of  a failure  in  the  forecast  of  the  large-scale 
circulation  pattern.  The  relative  success  of  the  operational  forecasts  in  these 
respects  may  have  been  due  to  a subjective  anticipation  of  this  behavior  on 
tin'  part  of  the  hurricane  forecasters  who  provided  the  operational  initial 
truck  lirections  and  speeds.  We  do  not  believe  that  this  result  would  be  found 
in  a sample  of  forecasts  in  which  more  of  the  actual  tracks  were  irregular. 

1 1 )F  NTIFIC'ATION  OF  THE  PREDICTED  STORM  CENTER 

The  nominal  position  of  Ihe  storm  center  in  the  SANRAR  forecast  has 
been  taken  as  the  mean  of  the  positions  of  the  absolute  vorticity  maximum 
and  the  stream-function  minimum,  when  both  were  present.  If  the  latter 
•lanishes,  as  it  often  does  due  to  truncation  error  or  strengthening  of  the 
large-scale  ilow  in  the  vicinity  of  the  storm,  then  the  position  of  the  storm 
identified  with  that  of  the  vorticity  maximum.  Tn  this  case,  the  forecast 
' ark  displays  a spurious  jog  to  the  right.  Rarely,  the  forecast  vorticity 
1 For  that  matter,  the  statistical  significance  of  the  improvement  in  Table 
4 is  not  large,  but  we  choose  to  be  guided  by  physical  reasoning  and  believe 
ilm  result.  We  would  bet  on  the  next  sample. 
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maximum  cannot  be  unambigously  identified;  then  the  storm  is  presumed  to 
have  lost  its  identity.  This  did  not  happen  in  the  1975  sample,  the  predicted 
storms  being  more  persistent  than  the  observed  ones,  especially  after  land- 
fall. In  all  cases  the  predicted  storm  position,  however  obtained,  was  ad- 
justed by  a vector  correction  appropriate  to  the  vector  discrepancy  between 
the  specified  initial  position  and  the  position  which  emerged  from  the  anal- 
ysis and  initialization  procedure. 

It  occurred  to  us  that  since  the  predicted  position  of  the  stream-function 
minimum  may  be  strongly  affected  by  truncation  error,  and  since  the  physi- 
cal basis  of  the  forecast  is  conservation  of  absolute  vorticity,  the  position  of 
the  storm  should  be  identified  with  that  of  the  absolute  vorticity  maximum  at 
all  times.  Strictly  speaking,  the  relative  vorticity  maximum  is  the  relevant 
entity  and  the  absolute  maximum  will  be  found  slightly  to  the  north.  The 
strength  of  the  earth-vorticity  field,  however,  is  utterly  small  compared  to 
that  of  the  relative-vorticity  field,  so  that  the  discrepancy  is  undoubtedly 
negligible.  Moreover,  the  latter  is  not  explicitly  carrier  in  the  forecast. 
Therefore,  we  are  satisfied  with  the  position  of  the  absolute  maximum. 

Accordingly,  for  the  rerun  predictions,  the  position  errors  were  recom- 
puted on  the  basis  of  the  vorticity  maximum  alone,  with  results  shown  in 
Table  7.  The  adjustment  in  this  case  was  based  on  the  small  discrepancy 
between  the  specified  storm  position  and  the  initialized  vorticity  maximum. 
Surprisingly,  the  errors  are  slightly  larger,  overall  and  for  most  storms 
and  ranges.  Since  the  vorticity  maximum  will  lie  to  the  right  of  the  nominal 
position,  we  might  expect  that  the  left  bias  in  the  original  rerun  forecasts 
would  be  reduced.  The  mean  algebraic  errors  of  the  two  sets  of  rerun  fore- 
casts can  be  compared  from  the  data  in  Table  8.  We  find  that  indeed  the 
left  bias  is  reduced  in  the  first  24-hour  period,  overall  and  for  five  of  the 
seven  storms.  For  the  second  period,  however,  the  overall  left  bias  is  in- 
creased, a result  due  toEloise  and  Gladys.  The  slow  bias  in  speeds,  more- 
over, is  slightly  increased. 

We  believe  that  these  curious  results  were  produced  by  cancellation  of 
errors  in  the  rerun  forecasts  based  on  the  nominal  position,  and  that  the 
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TABLE  7.  Mean  position  errors  (nautical  miles) 
for  SANBAR  reruns  based  on  vortieity  maximum 
alone . 


00  hr 

12  hr 

24  hr 

36  hr 

48  hr 

72  hr 

Amy 

0 

65 

108 

154 

245 

422 

(7) 

(6) 

(5) 

(5) 

(5) 

(3) 

Blanche 

0 

31 

86 

122 

_ 

(4) 

(3) 

(2) 

(1) 

Carol  ine 

0 

24 

39 

65 

198 

(5) 

(4) 

(3) 

(2) 

(1) 

Doris 

0 

43 

78 

117 

174 

305 

(10) 

(9) 

(8) 

(7) 

(6) 

(4) 

E loise 

0 

42 

73 

109 

185 

320 

(15) 

(14) 

(13) 

(12) 

(ID 

(9) 

Faye 

0 

64 

141 

182 

202 

421 

(14) 

(13) 

(ID 

(9) 

(7) 

(5) 

Gladys 

0 

56 

125 

203 

296 

437 

(17) 

(17) 

(16) 

(15) 

(14) 

(12) 

Bailie 

0 

100 

(2) 

(1) 

Mean 

0 

51 

103 

154 

229 

385 

(74) 

(671 

(58) 

(51) 

(44) 

(33) 

Change  from 

original 

reruns 

0 

+ 1 

+4 

+2 

+5 

+9 

( ) Number  in  sample 
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biases  in  these  forecasts  would  have  been  larger,  and  the  position  errors 
also  larger-,  than  in  forecasts  based  on  the  vorticity  maximum  alone,  had 
not  this  cancellation  occurred  . How  this  comes  about  is  shown  in  Fig.  3. 

The  initial  correction  for  the  nominal -position  forecasts  is  toward  the  north- 
east, on  the  average,  because  the  average  storm  is  embedded  in  a southeast- 
erly large-scale  flow  and  the  initialized  stream-function  minimum  is  too  far 
to  the  southwest.  When  this  minimum  disappears  in  the  course  of  these  fore- 
casts, as  more  than  half  did  by  48  hours,  the  position  (now  based  solely  on 
the  vorticity  maximum)  receives  an  undeserved  boost  toward  the  northeast, 
tending  to  reduce  the  average  bias  in  the  forecast  position.  The  vorticity - 
maximum  forecasts  do  not  receive  this  benefit  because  the  initial  discrep- 
ancy, and  thus  the  correction,  is  small.  When  the  differences  in  the  corr- 
ections for  the  nominal-position  and  vortieity-maximum  forecasts  (seen  in 
Pig.  3)  are  compared  with  the  differences  shown  in  Table  7,  it  appears  that 
the  success  of  the  nominal-position  forecast';  was  probably  a happy  procedural 
accident. 

Needless  to  say,  this  result  could  not  be  expected  in  a sample  of  fore- 
casts containing  relatively  fewer  recurving,  accelerating  tracks;  and  we 
recommend  use  of  the  vorticity  maximum  alone,  because  it  makes  better 
physical  sense.  Finally,  however,  the  entire  matter  is  not  very  important; 
the  ten-mile  reduction  of  error,  which  is  the  most  that  could  be  expected, 
would  have  little  impact  on  the  state  of  the  art.  The  major  areas  of  pros- 
pective improvement  lie  elsewhere. 

SANBAR  ANALYSIS  IN  OCEANIC  REGIONS 

Away  from  regions  of  adequate  rawinSonde  coverage,  the  analysis  of 
the  initial  large-scale  field  of  motion  is  not  only  difficult,  but  also  crucially 
important  for  most  tropical-storm  forecasts.  In  the  SANBAR  grid  area  an 
array  of  44  "hogu*  points”  was  established  (Fig.  4)  by  subjective  appraisal 
of  need.  Estimates  of  the  mean  tropospheric  wind  at  these  points,  however 
obtained  , were  regarded  as  having  the  same  statistical  properties  as  rawin- 
sonde  wind  observations,  and  were  therefore  used  in  the  construction  of 


Initialized  Vorficity  Maximum 


Specified  Position  J 


A ^ Correction  for  "Vorticfy  Maximum" 
Forecasts 


/ 5 


/ 

Initialized  Nominal  Position  J 


Correction  for  "Nominal  Posilion' 
Forecasts 


Initialized  Stream-Function 
Minimum 


24  hr 


+ 


48  hr 


+ 


72  hr 


l isurr  3^  Sketches  of  mrrpctinn  for  "nominal  position"  forecasts 
am!  for  "vorticify  maximum"  forecasts,  with  actual  values  of  mean 
discrepancies  given  in  nautical  miles  See  text. 
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regression  equations  for  performing  the  analysis  at  the  relevant  points  of  the 
SANBAR  computational  grid. 

Initially  ( in  1968),  the  bogus-point  values  were  obtained  subjectively  by 
a judicious  blend  of  the  current  operational  surface  and  200-mb  analyses  pre- 
pared manually  at  N11C  and  the  500-mb  analysis  (at  higher  latitudes)  prepared 
by  NMC,  occasionally  seasoned  with  winds  observed  by  reconnaissance  air- 
craft. These  information  sources,  save  the  last,  were  generally  12  hours 
old  at  the  time  the  bogus-point  estimates  were  made.  The  procedure,  more- 
over, was  awkward  and  time-consuming. 

In  time,  the  process  of  analysis  of  the  flow  pattern  at  low  levels  (ATOLL) 
and  near  200-mb  was  automated  (Wise  and  Simpson,  1971),  at  NHC;  the  auto- 
mation of  bogus-point  estimation  was  soon  to  follow  ( Pike,  1972).  Further 
and  most  importantly,  a new  data  source  became  available:  estimates  of 
wind  based  on  cloud-motion  vectors  obtained  from  successive  views  from 
geosynchronous  satellites.  The  cloud  elements  were  morphologically  char- 
acterized as  "low",  "middle",  or  "high".  The  first  and  last  were  added  to 
the  data  base  for  the  ATOLL  and  200-mb  analyses,  respectively.  The  middle- 
cloud  motions,  which  were  rarely  obtained,  were  not  used  systematically 
at  Nil C. 

Pike  (1975)  derived  regression  equations  for  estimating  the  tropospheric 
mean  wind,  required  by  SANBAR,  using  wind  observations  at  850  mb  "nd 
at  200mb  as  predictors.  In  the  small  dependent  data  sample,  rawinsonde 
observations  at  these  levels  were  presumed  to  simulate  winds  derived  from 
ATOLL  and  200-mb  analyses  in  actual  application.  Adams  and  Sanders!  1975) 
established  comparable  regression  equations  from  a very  large  sample  of 
data,  from  the  Pacific  as  well  as  the  Atlantic  sector,  stratified  by  geogra- 
phical area  and  by  season.  Their  results,  for  the  area  corresponding  most 
closely  to  Pike's  are  shown  in  Table  9,  along  with  his  earlier  results.  As 
can  be  seen,  they  found  little  advantage  in  stratifying  by  season. 

In  view  of  the  difference  between  these  definitive  equations  and  Pike's  , 
and  in  view  of  lack  of  knowledge  of  the  adequacy  of  the  automated  ATOLL 
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TABLE  9.  Regression  equations  for  zonal  and 
meridional  components  (knots)  of  mean  tropospheric 
wind,  based  on  rawinsonde  observations  at  850  mb 
and  200  mb. 

Pike  (June  - November) 


100  = -0.512  + 0.561  u850  + 0.399  o200 


100  ■ 0574  + °-269  ''850  + °'265  V200 


Adams  and  Sanders!  June -October ) 


100  ""  0,394  + °-530  u850  + °’374  u250 


■100  = -°-513  + °-450  v850  + 327  V250 


Adams  and  Sanders  (June -August) 


= 0.268  + 0.538  u_  + 0.  355  uocri 
100  850  250 


100  = -0.515  + 0. 437  v + 0.  322  vOK- 

850  250 


Adams  and  Sanders  (September-October) 


,00  = 0.581  - 0 521  u850  + 0.386  u250 


100  = -0.451  + 0.467  v850  + 0.332v,50 


n 


* 


E 


and  200-mb  analyses,  we  decided  to  reanalyze  the  oceanic  data  base  for  a 
selection  of  cases  from  the  1975  season.  The  original  intent  was  to  obtain 
revised  values  at  the  bogus  points  at  these  two  levels,  and  to  use  Adams's 
and  Sanders's  equations  to  estimate  the  tropospheric  mean.  The  list  of  cases 
is  given  in  Table  10.  We  initially  selected  18  bad  forecasts  and  18  good  ones 
from  the  operational  sample,  expecting  that  the  bad  ones  would  show  improve- 
ment and  the  good  ones  deterioration;  and  hoping  that  the  improvement  would 
exceed  the  deterioration.  The  first  result  would  likely  ensue  even  if  our  alter- 
ations had  no  real  merit,  while  the  second  would  assure  that  we  had  in  fact 
done  something  helpful.  In  the  event,  one  of  the  good  forecasts  could  not  be 
re run. 

Upon  our  realization,  however,  that  the  ATOLL  and  200-mb  analyses 

are  now  performed  exactly  on  the  SANBAR  computational  grid,  the  use  of 

2 

bogus  points  seemed  to  make  little  sense  , Unless  the  NHC  analysis  algo- 
rithms were  egregiously  bad,  more  useful  information  could  be  brought  to 
bear  by  applying  the  regression  equations  directly  to  the  analyzed  data  at 
the  SANBAR  computational  grid.  Accordingly,  we  are  carrying  out  the  re- 
vised forecasts  on  the  basis  of  a revised  analysis  technique  which  makes  no 
reference  to  the  bogus  points  and  uses  data  directly  from  the  ATOLL  and 
200-mb  analyses  in  regions  sufficiently  removed  from  rawinsonde  observa- 
tions. The  results  of  these  revisions  will  be  discussed  in  a later  scientific 
report. 


L We  feel  that  the  distinction  between  Pike's  use  of  200-mb  winds  and 
Adams's  and  Sanders's  use  of  250  mb  is  of  little  consequence.  The  NHC 
200-mb  analysis,  based  primarily  on  aircraft  observations  over  the  oceans 
probably  more  nearly  represents  250  mb  anyway. 

2.  Except  poleward  of  45°N,  the  northern  limit  of  the  ATOLL  and  200- nib 
analyses . 
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Radius  of  influence  assumed  300nm  in  all  cases  except  Gladys,  09^28  , 1200GMT,  and  Gladys,  09 '29.  0000GMT. 
for  which  radius  assumed  200nm. 
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DIAGNOSIS  OF  ERROR  IN  SOME  1975  FORECASTS 

Prior  to  analyzing  the  revised  forecasts  described  above,  we  undertook 
a study  of  the  particularly  poor  operational  forecasts  listed  in  Table  10,  with 
the  aim  of  delermining  the  causes  of  large  forecast  error  and  of  estimating 
what  improvement  might  be  expected  from  a better  analysis  of  the  initial 
large-scale  flow  pattern.  We  found  a variety  of  circumstances  leading  to 
especially  bad  forecasts,  the  major  ones  occurring  in  the  case  histories  to 
be  discussed  below. 

Amy,  July  3,  0000GMT 

This  storm  (see  Fig.  1)  originated  as  a depression  over  the  Gulf  Stream 
east  of  Florida  and  reached  tropical  storm  strength  east  of  Cape  Hatteras  on 
July  1.  It  reached  hurricane  strength  only  briefly,  if  at  all,  before  acceler- 
ating northeastward,  losing  its  identity  east  of  Newfoundland  on  the  4th.  Amy's 
main  effect  was  to  harass  yachts  competing  in  the  eastbound  transatlantic  race. 
And  to  damage  yachts  cruising  between  Bermuda  and  the  United  States,  taken 
unaware  by  this  early-season  storm  of  semi-tropical  character.  (Hebert,  1976). 
(According  to  Van  Gemert  H977),  winds  about  40nm  north  of  the  center  late 
on  the  3rd  were  between  50  and  60  knots,  with  seas  to  25  or  30  feet.  Condi- 
tions were  evidently  somewhat  worse  closer  to  the  storm  center]. 

The  details  of  the  predicted  and  observed  tracks  appear  in  Fig.  5.  The 
operational  track  direction  and  speed  at  the  initial  time  was  toward  065°  at 
10  knots.  The  storm  was  predicted  to  move  in  approximately  this  direction 
and  to  accelerate  only  slightly;  14  knots  in  the  first  24  hours  and  17  knots  in 
the  second.  1 In  tl  e event,  the  storm  curved  slightly  to  the  left  and  acceler- 
ated dramatically  to  24  knots  12  to  24  hours  after  initial  time.  The  result 
was  a large  205-nm  error  in  the  forecast  position  at  24  hours  range,  with 


1.  The  jog  of  the  forecast  track  to  the  right  between  36  and  48  hours  is  an 
artificial  result  of  the  loss  of  stream-function  mmimum,  as  discussed 
earlier . 
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commensurate ly  large  errors  at  later  times  had  we  not  received  the  blessing 

I of  dissolution  of  the  actual  storm.  The  best-track  position  and  track  velocity 
were  close  to  the  operational  values  and  would  have  led  to  no  great  improve- 
ment in  this  forecast,  but  we  provided  an  erroneous  velocity  toward  075°  at 
18  knots  for  the  rerun,  yielding  the  modified  track  shown.  The  storm  failed 
to  accelerate  through  the  first  48  hours  of  the  forecast.  (This  case,  of  course, 
was  not  included  in  the  sample  of  74  discussed  in  earlier  sections.  ) The  fail- 
ure of  either  the  operational  or  the  spurious  rerun  forecasts  to  show  early 
acceleration  led  us  to  believe  that  the  initial  large-scale  analysis  may  have 
been  seriously  in  error. 

This  flow  pattern,  from  the  operational  SANBAR  run,  appears  in  Fig.  6. 
Filtering  of  the  storm  leaves  a large-scale  flow  toward  the  east-northeast  at 
about  10  knots  over  the  initial  position  of  Amy,  while  large-scale  confluence 
produces  a speed  of  about  20  knots  at  the  24-hour  forecast  position.  * Thus, 
the  first  24-hours  of  the  forecast  track  was  a direct  consequence  of  the  initial 
analysis;  its  temporal  evolution  was  of  little  consequence. 

The  data  bases  for  the  ATOLL  and  200-mb  analyses  (supplemented  by 
some  data  evidently  not  available  in  real  time)  are  shown  in  Fig.  7.  In  the 
crucial  area,  say,  bounded  by  the  center  of  storm,  bogus  points  7 and  5, 
Newfoundland  and  Nova  Scotia,  there  are  few  observations  to  guide  the 
ATOLL  analysis.  At  200mb  , a number  of  reports  of  strong  southwesterly 
winds  lie  in  this  area,  but  not  at  bogus  points  5 and  7. 

The  strong  wind  observed  by  the  ship  at  37°N  59°W  and  the  strong  south - 

iy  southwesterlies  at  200  mb  derived  from  satellite  observations  near  37°N 

* o 

57  w as  well  as  the  estimated  value  at  bogus  point  4 would  have  been  ignored 
since  they  lie  within  the  maximum  influence  distance  of  the  storm.  It  is 

* clear  from  the  storm  parameters  (Table  10)  that  the  synthetic  analysis  at 
37°N  59°W  would  yield  a much  weaker  mean  wind  than  any  reasonable  re - 

f- 

1.  A useful  rule  of  thumb  is  that  the  wind  speed  in  knots  is  10  times  the 

number  of  stream-function  channels  contained  in  a distance  of  300  nm 

( 5 degrees  of  latitude). 

II 
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Figure  6.  Initial  analysis  of  large  scale  flow  pattern.  July  3.  0000GMT 
Solid  lines  are  isopliths  of  stream  function,  at  intervals  of  300x10  m"  s" 
from  the  operational  SANBAR  analysis.  Dashed  lines  are  subjectively 
analyzed  500-mb  contours  at  intervals  of  6 dm.  Triangles  indicate  500- mb 
heights  obtained  from  satellite  data-,  circles,  from  rawinsonde  observation* 
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Figure  7.  Data  bases  for  oceanic  analysis,  July  3.  0000GMT 
a)  ATOLL;  b)  200-mb  Positions  of  bogus  points  are  shown 
by  circled  x's  Radius  of  influence  is  shown  by  the  dashed 
circle.  Dashed  arrows  indicate  data  not  in  NHC  collection. 
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gression  estimate  based  on  direct  analysis  of  the  observed  data.  Inspection 
of  Fig.  5 confirms  this  expectation.  At  42°N  55°W  a combination  of  the  wind 
observed  by  aircraft  and  an  estimate  of  140°  at  15  knots  estimated  from  the 
ATOLL  data  would  yield,  upon  application  of  Adams's  and  Sanders's  .June- 
October  regression  equations,  a tropospheric  mean  wind  of  207°  at  22  knots. 

An  analysis  based  only  on  rawinsonde  and  bogus  observations  would  not  de- 
tect this  wind;  and  Fig.  5 shows  a direction  of  about  26Q1  at  this  point  In 
all,  the  initial  large-scale  analysis  appears  to  be* seriously  in  error,  to  the 
considerable  detriment  of  the  forecast. 

In  this  case,  it  further  appears  that  constant-pressure  height  data,  stu- 

2 

diously  avoided  in  the  tropically-oriented  SANBAR  system,  might  have  been 
helpful.  Fig.  8 contains  the  heights  at  the  500-mb  level,  which  yields  a geo- 
strophic  wind  in  close  approximation  to  the  tropospheric  mean  wind  In  add- 
ition to  values  at  the  rawdnsonde  stations,  a number  of  estimates  were  made 
by  adding  a 1000-mb  he ight  derived  from  the  sea-level  pressure  analysis  to 
the  thickness  of  the  layer  from  lOOOmb  to  500mb  derived  from  SIRS  temper- 
ature soundings  from  the  NOAA-5  polar  orbiter.  A subjective  analysis  of 
the  500-mb  contours,  based  on  these  data,  appears  in  Fig.  8 ( a > ( Though  the 
satellite-derived  heights  are  clearly  poorer  than  the  rawinsonde  observations, 
the  time  series  of  rawinsonde  data  at  Bermuda  , near  32°N  65°W  . suggests 
that  the  observed  height  is  erroneously  low  at  this  time  and  that  the  satellite 
data  are  essentially  correct.  ) The  subjective  analysis  yields  geostrophic 
winds  which  are  in  accord  with  our  foregoing  estimates  based  on  observed 
wind  data.  A comparison  of  the  wind  fields  derived  from  the  subjectively 

1.  The  large  discrepancy  between  the  stream  flow  and  the  observed  wind  at 
Shelburne,  Nova  Scotia  (near  44(  N 65  W)  is  due  to  the  omission  of  this  re- 
cently-installed station  from  the  list  used  in  the  SANBAR  regression  equations. 
This  list  and  the  dependent  equations  should  be  updated  immediately. 

2.  This  is  a slight  overstatement,  since  the  ATOLL  and  200-mb  analyses 
use  NMC  grid-point  wind  values  ab  a first  guess,  and  these  have  been  derived 
in  part  from  rawinsonde  and  satellite  heights. 
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Figure  H Wind  fields  from  a)  the  subjective  500- mb  ger 
strophic  analysis  and  bl  subjectively-smoothed  SANBAF 
initial  analysis  Dotted  lines  indicate  isotachs  at  10-kt 
intervals  Other  notation  the  same  as  Figure  6 


smoothed  SANBAR  initial  analysis  and  from  the  500-mb  geostrophic  analysis 
appear  in  Fig.  8. 

We  observe,  incidentally,  from  the  ATOI.L  data  that  the  storm  influenced 
a region  with  a radius  of  about  500  nm  rather  than  the  nominal  300-nm  value 
used  in  the  SANBAR  calculation.  This  problem,  however,  is  interrelated 
with  the  yet  unresolved  problem  of  how  to  deal  effectively  with  data  within  the 
influence  region  of  the  storm,  as  discussed  by  Sanders  et  al  . (1975).  These 
important  problems  are  not  addressed  in  our  present  work. 

We  conclude  that  the  poor  forecast  of  the  track  of  Amy  was  due  in  this 
instance  to  faulty  initial  analysis  of  the  large-scale  flow  surrounding  the 
storm. 

Faye,  September  26,  0000GMT 

This  storm  was  first  detected  as  a typical  westward-moving  tropical  de- 
pression in  the  central  Atlantic  on  September  18.  The  track  in  Fig.  1 shows 
that  Faye  barely  reached  hurricane  strength  by  the  initial  time  of  this  fore- 
cast, after  a long  gestation  period.  She  had  begun  recurvature  on  the  24th  , 
moving  slowly  and  somewhat  erratically.  With  acceleration  and  recurvature 
into  the  middle-latitude  westerlies  far  from  land,  Faye  represented  only  a 
threat  to  shipping.  This  threat  may  have  been  considerable,  however,  since 
her  abrupt  eastward  acceleration  on  the  28th  brought  her  path  to  speeds  which 
only  the  fastest  and  strongest  ships  could  out-maneuver. 

The  SANBAR  and  actual  tracks  are  shown  in  detail  in  Fig.  9,  for  the  pre- 
diction starting  at  0000GMT  on  the  26th.  Position  errors  for  both  the  opera- 
tional and  the  best-track  rerun  forecasts  are  given  in  Table  11. 


TABLE  11.  Position  errors  (nautical  miles)  in  SANBAR 
forecasts  for  Faye,  September  26.0000GMT. 

Initial  track 


Forecast  Dir. 

(°)  and  Spd(kts) 

0 hrs 

24  hrs 

48  hrs 

72  h rs 

Operational 

315 

-11 

21 

212 

345 

961 

Best  track  rerun 

330 

-11 

0 

183 

253 

772 

Change  from  operational 

-21 

-29 

-92 

-189 

The  best-track  information  reduced  the  position  errors  considerably,  with 
the  improvement  increasing  monotonically  to  72  hours.  Reference  to  Fig.  9, 
however,  shows  that  the  24-hour  benefit  was  due  solely  to  the  improved  track 
direction  in  the  best-track  forecast,  and  that  the  reduced  errors  at  later  times 
were  attributable  mainly  to  the  preliminary  recurvature  in  the  rerun  forecast. 
Speed  was  a major  problem  in  both  cases:  both  the  immediate  jump  from  the 
initial  11  knots  to  an  actual  17  knots,  and  the  observed  mighty  increase  in 
speed  from  48  to  72  hours.  Thus,  we  believe  that  the  large-scale  flow  was 
both  incorrectly  analyzed  initially  and  unusually  badly  predicted. 

The  initial  stream-function  pattern  appears  in  Fig.  10,  while  the  data 
bases  f6r  the  ATOLL  and  200-mb  analyses  appear  in  Fig.  11.  At  bogus  points 
9,  13,  and  14,  at  Bermuda,  and  at  20°N  60°W  it  was  possible  to  compare 
large-scale  winds  derived  from  the  stream  pattern  with  winds  computed  from 
subjective  estimates  from  the  ATOLL  and  200-mb  data  bases.  These  esti- 
mates agreed  with  operational  values  at  some  but  not  all  of  these  points.  In 
cases  of  disagreement,  an  observation  was  close  to  the  point.  It  can  be  seen 
from  the  results  in  Fig.  10  that  the  directions  of  the  stream-function  and  es- 
timated winds  agree  closely  except  at  bogus  point  14,  where  the  latter  is 
substantially  veered  from  the  former.  The  estimated  speeds,  however,  are 
greater  without  exception,  substantially  so  near  and  perhaps  south  of  the  cen- 
ter. These  speeds,  in  fact,  correspond  in  the  mean  closely  to  the  average 
of  the  actual  track  speeds  six  hours  prior  and  12  hours  after  the  initial  time. 
(A  correct  initial  analysis,  after  all,  should  provide  a goodhindcast  as  well 
as  a good  forecast.  ) There  is  some  indication  of  a downstream  increase  in 
speed,  suggesting  acceleration  along  the  track,  but  not  the  sudden  speedup 
ob  se  r ve  d . 

The  question  then  arises  why  the  SANBAR  predicted  speeds  in  the  first 
24  hours  were  so  slow.  First,  winds  derived  from  observations  at  bogus 
points  9 and  13  were  discarded,  since  they  lie  within  the  outer  portions  of 
the  300-nm  radius  deemed  to  be  influenced  by  Faye.  They  were  replaced 
by  weaker  winds  based  on  the  11-knot  initial  track  speed.  It  appears,  then, 
that  the  present  SANBAR  analysis  algorithm  loses  valuable  information  in 
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Figure  11.  Databases  for  oceanic  analysis,  September 
26.  0000GMT;  a)  ATOLL:  b)  200-mb.  Notation  the 
same  as  Figure  7. 
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the  outer  portion  of  the  region  influenced  (only  slightly)  by  the  storm. 

Second,  it  can  be  seen  from  Fig.  11  that  the  large-scale  flow  displayed 
substantial,  well -organized,  shear  between  the  lower  and  upper  troposphere 
in  the  region  surrounding  Faye.  An  equivalent-barotropic  effect  could  easily 
be  added  to  the  SANBAR  forecast,  but  as  pointed  out  by  Sanders _et  _ad.  (197  5), 
it  is  in  general  neither  substantial  nor  reliable  in  the  tropical  regions  in  which 
this  model  is  currently  employed.  (This  effect,  if  applied  in  the  present  in- 
stance, would  yield  a better  short-term  forecast  but  a poorer  hindcast;  on 
another  occasion  it  might  be  the  other  way  around.  ) Fig.  11  also  shows  that 
Faye  was  embedded  in  a large  trough  system  in  the  ATOLL  layer.  This  pre- 
sumably same  trough  is  evident  in  the  200- mb  data,  with  stronger  circulation 
and  substantial  southward  displacement.  There  is  undoubtedly  positive  ther- 
mal-vorticity  advection  over  Faye,  similar  to  developing  extratropical  cyclones, 
however  modest  the  magnitude  in  this  case.  A fully  baroclinic  model  would 
be  required  to  calculate  accurately  the  consequence  of  this  structure. 
Hovermale's  MFM  Model  (1975)  should  be  helpful  in  this  respect. 

Third,  as  can  be  seen  from  Fig.  10,  there  are  no  satellite  data  to  the 
right  of  the  path  of  the  storm.  They  might  have  been  helpful,  as  in  the  case 
of  Amy  discussed  above. 

In  summary,  there  were  obvious  flaws  in  the  initial  analysis,  which  were 
damaging  to  the  short-term  forecast.  Part  of  the  answer  remains  shrouded 
in  mystery. 

The  failure  of  the  model  to  produce  the  dramatic  forward  acceleration 
observed  from  48  to  72  hours  is  another  story,  only  part  of  which  lies  in  the 
failure  of  SANBAR  to  put  Faye  in  the  right  place  at  48  hours.  The  matter  is 
elucidated  in  Fig.  12,  in  which  the  large-scale  predicted  and  observed  48-hour 
stream- function  changes  are  depicted.  Since  the  stream-function  field  can 
be  specified  only  as  the  residual  from  an  arbitrary  constant,  some  assump- 
tion must  be  made  to  arrive  at  these  fields.  We  assumed  that  the  stream- 
function  value  remained  essentially  unchanged  in  the  lower  right-hand  corner 
of  the  grid  area. 
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We  note  that  the  intensity  of  the  forecast  changes  is  not  as  great  as  that 
of  the  observed  ones.  This  failure  is  probably  typical  of  any  barotropic 
model.  Tn  particular,  the  predicted  gradient  of  change  in  the  path  of  Faye 
was  woefully  in  error.  The  observed  maximum  rise  along  25°N  in  the  At- 
lantic was  only  hinted  at  in  the  forecast.  There  were  substantial  changes 
along  the  boundaries,  held  approximately  constant  in  the  SANBAR  forecast, 
but  these  were  probably  not  primarily  to  blame,  since  the  centers  of  change 
were  well  removed  from  the  limits  of  the  forecast  grid.  Note  particularly  the 
prominent  center  of  fall  along  the  St.  Lawrence  River,  due  to  the  advance  of 
the  pronounced  trough  initially  in  the  Mississippi  Valley.  This  prominent 
change,  feebly  simulated  in  the  SANBAR  forecast,  was  doubtless  responsi- 
ble, along  with  the  general  rise  along  25°N  , for  the  observed  increase  in 
speed  of  the  storm.  The  large-scale  flow  speed  in  the  path  of  Faye  was  25- 
30  knots  and  increasing  with  time. 

The  large-scale  NMC  forecasts  for  approximately  this  time  appear  in 
Fig.  13*.  It  is  clear  that  neither  the  barotropic  ncr  the  six-layer  primitive- 
equation  (P'S)  forecast  adequately  predicted  the  substantial  changes  in  the 
large-scale  flow.  The  PE  72-hour  forecast  did  show  the  trend,  but  too  little 
and  too  late.  From  this  we  conclude  that,  as  suggested  above,  fixed  boundary 
values  of  stream  function  were  not  a major  contributor  to  the  SANBAR  fail- 
ure in  this  case,  and  that  although  baroclinic  effects  were  doubtless  impor- 
tant they  were  not  adequately  represented  in  the  PE  forecast. 

In  summary,  the  large  errors  at  48  and  72  hours  were  due  to  egregious 
error  in  the  large-scale  SANBAR  forecast,  principally  due  to  the  'ailure 
to  consider  baroclinic  effects  and  only  secondarily  to  the  failure  to  provide 
for  boundary  changes. 


sets  of  NMC  forecasts  were  not  available  at  MIT  at  this 


1.  Comple 
time. 
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Figure  13.  NMC  prognostic  charts:  a)  36-hour  barotropic 
forecast  valid  at  0000GMT,  September  27:  b)  24-hour  PE 
forecast  valid  at  OOOOGMT.  September  27:  and  c)  72-hour 
PE  forecast  valid  at  OOOOGMT.  September  2!' 
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Gladys,  September  30.  0000  GMT 

As  shown  in  Fig.  1,  Gladys,  like  Faye,  grew  from  a tropical  depression 
that  doubtless  evolved  from  an  easterly  wave  emerging  from  Africa  (Burpee, 
1972).  Gladys,  however,  achieved  hurricane  strength  much  more  quickly 
and  followed  a track  somewhat  south  and  west  of  Faye's,  although  generally 
parallel  to  it.  The  former  storm,  in  fact,  can  be  seen  in  Figs.  10-12,  at  a 
time  when  serious  SANBAR  forecast  errors  had  not  yet  developed.  As  in 
the  Faye  forecast  discussed  above,  the  primary  failure  is  inadequate  pre- 
diction of  northeastward  acceleration  following  recurvature,  but  the  cause 
in  this  case  is  different,  as  we  shall  see. 

The  details  of  the  predicted  and  observed  tracks  following  the  initial 
time  of  0000GMT,  September  30  appear  in  Fig.  14.  , while  the  ATOLL  and 
200-mb  data  bases  appear-  in  Fig.  15.  Position  errors  for  the  operational 
and  rerun  forecasts  are  given  in  Table  12.  The  24  and  48-hour  forecasts 
in  this  case  are  better  than  average,  the  latter  spectacularly  so.  There  is 


TABLE  12.  Position  errors  (nautical  miles)  in  SANBAR 
forecasts  for  Gladys,  September  30,  0000GMT. 

Initial  track 


Forecast 

Dir.  (°)  and  Spd(kts) 

0 hrs 

24  hrs 

48hrs 

72  hrs 

t 

Operational 

295 

14 

5 

50 

15 

661 

Best  track  rerun  295 

14 

0 

47 

15 

675 

l 

Change  from  operational 

-5 

-3 

0 

+ 14 

no  significant  difference  between  the  operational  and  rerun  forecasts,  as  the 
position  and  track  velocity  were  well  known  in  real  time.  The  paucity  of 
initial  data,  particularly  at  200  mb  east  of  Gladys  (see  Fig.  15),  demon  - 
strates  that  one  can  be  lucky  on  occasion.  The  observations  in  Fig.  15(b), 
however,  are  adequate  to  show  the  upper-level  trough  wrapped  around  the 
western  and  southern  periphery  of  the  storm,  extending  thence  eastward 
between  latitudes  1 5 'N  and  20°N.  This  feature,  undoubtedly  the  same  as  the 
one  seen  in  Fig.  11(b),  does  not  seem  in  this  instance  (contrary  to  the  Faye 
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station  the  same  as  Figtir 


Figure  15  Data  bases  for  oceanic  analysis,  September  30. 
0000CMT:  a)  ATOLL;  b)  200-mb.  Notation  the  same  as 
Figure  7. 
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case  discussed  above)  to  be  providing  positive  vorticity  advection  over  the 
storm  area,  in  which  the  flow  is  nearly  symmetric.  The  difference  may  be 
a reflection  of  Faye's  growth  vs  Gladys's  maturity  at  the  time  studied. 

The  large  error  is  in  the  range  of  48  to  72  hours,  when  Fig.  14  shows 
that  the  predicted  storm  came  nearly  to  a halt  while  the  living  Gladys  leapt 
ahead  to  a speed  of  40  knots  in  the  last  half  of  the  interval.  The  flow  pattern 
and  48-hour  observed  changes  shown  if  Fig.  16(a)  display  a series  of  mobile 
ridges  and  trough  in  the  westerlies  poleward  of  35°N  with  considerable  am- 
plitude but  little  indication  of  change  of  intensity.  The  predicted  changes  in 
Fig.  16(b)  , on  the  other  hand,  show  little  of  this  mobility  but  rather  a sub- 
stantial average  stream-function  rise.  The  band  of  maximum  rise  extending 
north  of  the  storm  path  between  latitudes  40°N  and  45°N,  together  with  the 
westward  displacement  of  the  subtropical  anticyclone  from  an  initial  position 
at  31°N  56°W  to  a position  near  33°N  65°W  48  hours  later  explain  Faye's 
subsequent  predicted  halt'!  This  large-scale  anticyclone  doubtlessly  had 
settled  directly  over  the  storm  by  72  hours.  The  actual  storm,  on  the  other 
hand,  was  presumably  caught  up  in  the  strong  southwesterly  flow  between  the 
trough  in  the  Great  Lakes  region  and  the  downstream  ridge  over  the  Canadian 
Maritime  Provinces  at  48  hours  range. 

The  failure  of  the  SANBAR  forecast  in  this  case  has  two  causes.  First, 
the  observed  changes  were  large  along  the  fixed  northern  boundary  while  the 
maximum  magnitudes  were  not  far  from  it.  Second,  the  average  flux  of  earth 
vorticity  into  the  forecast  area  appears  to  be  in  the  anticyclonic  sense,  pri- 
marily because  of  the  influx  of  low  values  at  the  left  boundary  near  latitude 
25°N  and  the  efflux  of  large  values  at  the  right  boundary  near  40°N.  While 
the  forecast  program  is  being  altered  to  allow  for  stream-function  changes 
on  the  boundary,  it  might  be  a good  idea  to  allow  for  vorticity  changes  as 
well  . Of  course,  only  the  local  changes  of  relative  vorticity  would  count; 
their  importance  has  not  yet  been  demonstrated  but  is  probably  considerable 


I.  Only  the  initial  and  48-hour  predicted  .SANBAR  flow  patterns  are  arch- 
ived by  NHC. 
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on  occasion. 

Our  diagnosis  is  supported  by  the  barotropic  and  PE  forecasts  shown  in 
Fig.  17.  The  general  trend  of  events  is  successfully  indicated,  especially 
by  the  latle  r,  enough  so  to  produce  substantial  improvement  had  they  been 
used  to  predict  the  storm  path  but  probably  not  accuately  enough  to  yield  a 
really  good  72 -hour  forecast. 

CONCLUSIONS 

Our  study  of  SANBAR  forecasts  made  operationally  during  the  197  5 
hurricane  season  yields  the  following  conclusions  and  recommendations: 

1.  A more  accurate  initial  position  and  track  will  produce  a substan- 
tial increase  in  accuracy  of  the  position  forecasts  at  ranges  of  24 
hours  and  probably  at  48  hours  . 

2.  The  use  of  satellite-derived  heights  may  provide  useful  informa- 
tion for  prediction  of  storms  initially  poleward  of,  say,  30°N. 

3.  The  current  method  of  analyzing  the  initial  stream-function  field 
in  the  a"ea  influenced  by  the  storm  circulation  is  resulting  in  loss 
of  useful  observational  information  (mostly  satellite-derived  wind 
estimates)  in  the  outer  fringes.  This  method  should  be  altered. 

4.  The  present  use  of  SANBAR  bogus  points  is  an  unnecessarily  com- 
plex procedure  which  is  probably  losing  valuable  information  con- 
tained in  increasingly  abundant  satellite  information.  The  anal- 
ysis method  should  be  altered,  at  least  over,  oceanic  regions. 

5.  Fixed  values  of  stream-function  and  absolute  vorticity  on  the 
boundaries  of  the  current  SANBAR  forecast  area  are  producing 
serious  errors  in  some  forecasts  beyond  24  hours  range.  Either 
the  forecast  area  should  be  enlarged  or  predicted  values  from  a 
hemispheric  or  global  model  should  be  applied  to  the  SANBAR 
boundaries. 

6.  Some  storm  tracks  respond  to  baroclinic  effects  in  the  large- 
scale  features  of  the  surrounding  region.  In  some  instances  the 
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storm  track  may  be  influenced  by  asymmetric  baroclinic  effects 
in  the  storm  structure  itself.  A baroclinic  model,  appropriate  to 
the  character  of  the  oceanic  data  base,  should  be  developed. 

7.  The  McIDAS  system  should  be  helpful  with  respect  to  items  1,  2, 

3,  and  4. 

8.  Even  when  all  this  is  done,  we  shall  still  see  in  part  and  prophesy 
in  part.  But  we  may  hope  to  see  through  the  glass  less  darkly. 
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